Nonmonotonic behavior of resistance in a superconductor-Luttinger liquid junction 
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Transport through a superconductor-Luttinger liquid junction is considered. When the inter- 
action in the Luttinger liquid is repulsive, the resistance of the junction with a sufficiently clean 
interface shows nonmonotonic temperature- or voltage-dependence due to the competition between 
the superconductivity and the repulsive interaction. The result is discussed in connection with 
recent experiments on single-wall carbon nanotubes in contact with superconducting leads. 
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PACS numbers: 73.20.Dx, 73.40.Hm 

In one dimensional systems, it is believed Q that the 
usual Fermi liquid description breaks down and systems 
become a Luttinger liquid (LL). LL behavior has been 
predicted || for single- wall carbon nanotubes (CNT), 
and demonstrated in experiments (3). Recently hybrid 
systems, which consist of LLs and superconductors (S), 
have also received theoretical attention ||, ||, (?J |). 
A couple of recent experiments |)| [l(J reported transport 
measurements on CNTs in contact with superconductors. 

The contact with superconductors may generate inter- 
esting phenomena. In a S-normal metal (N) junction, 
the junction resistance may be reduced due to the prox- 
imity effect mediated by the Andreev process fll| . This 
effect is more evident in a junction with a clean interface. 
In a S-LL junction, an interesting competition may arise 
between the superconductivity and the electron-electron 
interaction since the repulsive interaction in a LL tends 
to enhance the resistance |Lq| . The competition was in- 
vestigated in various theoretical works. However previ- 
ous studies were mostly focused either on the tunneling 
regime Q, where the Andreev process is strongly sup- 
pressed, or on the low energy regime g, §, § far below 
the superconducting gap A, where the repulsive interac- 
tion always wins over the superconductivity. 

In this paper, we investigate the competition in the 
energy range ~ A for a S-LL junction with a clean inter- 
face. We first examine two extreme energy regimes, far 
below A and far above A, and interpolate these results to 
the intermediate energy regime. Although quantitative 
predictions cannot be obtained with this approach, we 
still find that qualitative behaviors at the energy scales 
~ A can be determined without ambiguity. In particular, 
it is found that the competition between the supercon- 
ductivity and the interaction results in a nonmonotonic 
behavior of resistance (Figs, [j] and ||) when the junction 
interface is sufficiently clean. The result is discussed in 
connection with a recent experiment JuJ. 

Model: A purely one-dimensional model of a S-LL junc- 



tion is used, where the superconductor is also modeled 
as a one-dimensional conductor with the superconducting 
potential A (]]]. We choose the S to reside in the nega- 
tive semi-infinite space x < and the LL in the positive 
semi- infinite space x > 0. The bosonization formula JjJ, 

^ ri<T ocexp irkFX+ i(r9 c + ro-9 s + (t) c + o-(j)s)/V2 , (1) 

is used, where r = ±1 for R (right)/L (left) and a = ±1 
for spin up/down. The bosonic fields satisfy the commu- 
tation relation [0n(x),d y (f>n'(y)] = —iir5(x—y)5 fJ _ fl i, where 
/j,, fj/ = c (charge) or s (spin). We study the bosonized 
Hamiltonian H — Hll + H$ + H imp . Hll is the usual 
LL Hamiltonian, 
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where K c < 1 (repulsive interaction) and K s = 1 (spin 
SU(2) symmetry). Hs is given by Hn + Ha, where 



f° dx 

Ha = —A / — cos cos V2d s , 

J-oo x o 



(3) 



is the bosonized expression for J dx^ai^^R a^L -a + 
h.c), and xq is the short length scale cut-off. A = A(T) 
is assumed to be a real positive constant. Here ff/v is 
identical to Hll, except that K CjS = 1 and the integra- 
tion over x runs from — oo to 0. Finally the impurity 
backscattering potential at the interface is 



Hu 



-U cosV29 c (x = Q)coaV29 s (x = 0). (4) 



A similar Id model is suggested previously in Refs.[]|, 

Ideal clean interface: We first consider an ideal inter- 
face with Uq = 0. Above the critical temperature T c , 
where A(T) = 0, Ha vanishes and the system reduces 
to the N-LL junction with an ideal interface. By using 
the technique in Ref. p2[, one can integrate out degrees 
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grees of freedom for x 7^ and obtain an imaginary time 
effective action, 



FIG. 1: A schematic diagram of the zero bias resistance R as 
a function of the temperature T. The curve A shows R(T) for 
an ideal interface (Uo = 0). The curves B, C, and D describe 
progressively less ideal interfaces, and the curve E shows R(T) 
for a poorly transmitting interface (large Uo)- In the curves B 
and C, R(T) below T c is nonmonotonic due to the competition 
between the superconductivity and the repulsive interaction. 
The black squares in the curves B and C denote the crossover 
temperature T cr . J?nl = G^ and Rsl = Gg^ ( see text). 



of freedom for all x ^= and obtain an imaginary time 
effective action S for 9°(t) = 9 fJl (x = 0, r), 



S 



^ 7T/3 



K: 



|02(w„)| a + |0>n)| S 



(5) 



where (3 = T 1 is the inverse temperature, and 
0°(w) = f£ dTe^ T 0°{T). From the correlation func- 
tion (0°(t)0°(O)}, one obtains the zero bias conductance 
G = G NL = 2K C /(1 + K C )G @ where G = 2e 2 //i. 

On the other hand, for T < T c , where A(T) > T, H A 
becomes important and the cosine functions in Ha can be 
expanded around a potential minimum, which results in 
mass terms for the fields (j> c and 9 S . Since fluctuations of 
(f> c and 9 C are strongly suppressed by the mass terms for 
x < 0, and 6° may be regarded as constants, for exam- 
ple, zero. In the fermionic description, <f>® = 9 a s = sup- 
plies the boundary conditions (BC), ^rt(O) = — iip' L ,(0), 

ijjRl(0) = +iijj L *(0), which, combined with Hll, define a 
boundary conformal field theory in the semi-infinite space 
x > @. As demonstrated in Ref. fl, the BCs imply 
that all particles incident on the junction interface from 
the LL side are Andreev reflected, resulting in a trans- 
port of spin singlet Cooper pairs through the junction 
interface. 

The perfect Andreev reflection is expected to enhance 
the conductance. To evaluate G, we integrate out de- 
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9 C °K)| 2 +Q( W „)|0°K)| ; 



. (6) 



where Q(co) oc A/|w| for \u\ <C A. Fluctuations of 0" 
are massive as expected. On the other hand, fluctua- 
tions of 6® not only remain massless but are enhanced 
since the coefficient of |0°| 2 term is reduced compared to 
that in Eq. (ph. The enhancement is due to the An- 
dreev process and leads to the increased conductance 
G = Gsl = 2K c Go |l3[ l, which is larger than Gnl by 
the factor 1 + K c [|l3| , a generalization of the factor 2 en- 
hancement in a S-N junction with an ideal interface |jTl| . 
We then interpolate Gnl and Gsl to obtain the temper- 
ature dependence of R = G _1 (curve A in Fig. 

Almost ideal interface: When Uq is not zero but suffi- 
ciently small, the effective actions (||) and (||) are weakly 
perturbed by drH imp , which suppresses fluctuations 
of 0° and 0°, and reduces G. Above T c , this gener- 
ates a correction SG oc — T ai to Gnl, where a.\ = 
(K c — l)/(K c + 1) < 0, and far below T c , a correc- 
tion SG oc -T Q2 to G SL , where a 2 = 2(K C - 1) < 0. 
The power-laws can be verified by using the method 
in Ref. jl2j and the latter exponent is previously ob- 
tained in Refs. [?]]. The T-dependence of R can 
be constructed from the interpolation (the curve B in 
Fig. g). Note that while dR/dT < for T > T c and 
T <C T c due to SG, the interpolation requires dR/dT 
to be positive in the intermediate range T < T c since 
Gsl, Gnl ^ SG for a sufficiently clean interface and 
Gsl > Gnl- This result can be confirmed in an alter- 
native way. Direct calculation of G near T c , where both 
ifimp and Ha can be treated perturbatively, results in 
G = Gnl + c imp C/ 2 r Q "»p + c A A 2 T Q *, where a ifR = ati 
and the constant Cj mp < while a a — ~2 |L4] and 
ca > 0. Thus for a sufficiently clean interface (Uq — > 0), 
the temperature dependence of G is governed by the 
last term and dR/dT > since c A > §. This In- 
dependence of R for T < T c has not been revealed in 
earlier works and constitutes one of main results of this 
paper. 

For later use, here we define a crossover temperature 
T CI (< T c ), above which the superconductivity is domi- 
nant and below which the interaction effect is dominant, 
as the temperature where dR(T)/dT = (Fig. |l|). We 
also remark that at very low temperatures T -C T C r, the 
perturbative results near this "ideal-junction" fixed point 
(FP) are not valid since SG diverges at zero temperature. 
This regime is discussed below. 

Poorly transmitting interface: When the bare Uq = 00, 
the junction is disconnected and G = 0. Although the 
large bare Uq limit is not our concern, the diverging 
correction SG for a clean interface suggests that this 
"disconnected-junction" FP is a natural candidate for the 
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new FP that determines the resistance of a clean junction 
for T <C T cr . At this FP, we may impose the fermionic 
BC, ipL,(?(x = 0+) = e~ lS %jJn i<r (x = 0+), where S is a real 
constant, since all electrons incident on the junction in- 
terface should be normal-reflected. Combined with Hll, 
the BC defines another boundary conformal field theory 
in x > @. 

When the effective Uq is large but finite, current flow 
between the S and the LL can be generated by tunneling 
processes, which can be regarded as perturbations to 
the disconnected-junction FP. For T <C T CI (< T c ), the 
single-particle tunneling is exponentially suppressed. 
We thus consider only the Cooper pair tunneling H t = 



dV/dl 



#La(Q-)V'L-a(0-)^,-<r(0+)^fl ) a(0+) + h.C. 



which, upon bosonization, becomes 
H t = t cos V2« Li - < s ) cos \/2(0° 



(7) 



"H.LL/8 = Y>0 = 0± ) and e l,LL/S = 9 A X = 0± )' 



where 

One then obtains G oc T Q3 , where a 3 = 2(^ f T 1 - f ) > 0. 
Note that the power-law dependence justifies the ne- 
glect of the exponentially suppressed single-particle con- 
tribution. The same exponent is previously obtained in 
Refs. H M. All curves in Fig. |l] for nonzero bare Uo follow 
the power-law R oc T~ Q3 at sufficiently low temperatures, 
where the interaction wins over the superconductivity. 

For completeness, we also consider the temperature de- 
pendence above T c in a junction with sufficiently large 
bare Uo- Above T c , A = and the single-particle tun- 
neling leads to G oc T" 4 where a 4 = (R- 1 - l)/2 > 0. 
Note that and «4 are both positive, implying that the 
disconnected-junction FP is a stable one. The curve E in 
Fig. |l] schematically shows the resistance of the system 
when the bare Uo is sufficiently large. Note that R rises 
more rapidly below T c since a% > 0:4. 

General interface: To obtain qualitative behavior for 
arbitrary bare Uo, we interpolate the two opposite lim- 
its (curves B and E in Fig. |f|). The interpolation as- 
sumes the absence of a "partially-transmitting-junction" 
FP. The assumption is supported by a recent work 0, 
and also justifies the interpolation to the intermediate 
temperature range used to construct the curves in Fig. [|. 
A natural interpolation to intermediate Uq is depicted 
in Fig. [j]: For nonzero Uq larger than a critical value 
t/0,0 R increases monotonically with lowering tempera- 
ture (curves D and E) while for Uq < Uo, c , R shows non- 
monotonic behavior, that is, dR/dT > for T cr < T < T c 
and dR/dT < for T < T cr (curves B and C). 

dV/dl-V characteristic: Here we study the dV/dl- 
V characteristic at T < T c for almost ideal interfaces. 
We also confine ourselves to the subgap range eV < A. 
When T is slightly below T c , effects of the superconductor 
are weak and the characteristic should be similar to the 
characteristic for a N-LL junction, for which dV/dl de- 
cays monotonically with increasing V (curve 10 in Fig. |^). 
As T is lowered towards T cr , effects of the superconduc- 



10 
































3 



dV/dl 



A 

Fig. 2 Lee et al 



eV 



FIG. 2: A schematic drawing of the dV/dl-V characteristic 
for an almost ideal interface. T w T c for the curve 10 and the 
temperatures for the curves 9, 8, . . . , 1 are progressively lower. 
The curve denotes the characteristic at T = 0. Curves are 
vertically shifted for the clarity. Without the shift, all curves 
should almost collapse over each other near eV = A. When 
T is below the crossover temperature T CI (temperature for 
the curve 5), the characteristic becomes nonmonotonic in V 
due to the competition between the superconductivity and 
the repulsive interaction. The black squares in the curves 0-5 
denote the crossover voltage V C r(T). 



tor become stronger. Near the zero bias, according to 
the curves B and C in Fig. [I], the differential resistance 
should decay due to the effects of the superconductor. At 
relatively high bias eV » A, on the other hand, effects 
of the superconductor should remain still weak since the 
small energy barrier A — eV for single particle transport 
can be overcome by thermal fluctuations. As a result of 
different effectiveness of A at low and high biases, dV/dl 
develops a dip at zero bias (curves 9 to 5 in Fig. ||). This 
dip becomes deeper as T is lowered. 

When T is lowered below T cr , the large effective Uq 
due to the repulsive interaction wins over the effects of 
the superconductivity as shown in the curves B and C in 
Fig. |. Thus near the zero bias, dV/dl should grow as 
T approaches zero. However for sufficiently large bias, 
only minor changes in dV/dl are expected since the bias 
voltage itself serves as an energy cutoff for the RG flow of 
Uq and prohibits the crossover. As a result, in the tem- 
perature range < T < T cr , the characteristic becomes 
nonmonotonic in V: dV/dl has a negative slope below a 
certain temperature-dependent crossover bias V cr (T) and 
a positive slope above V cr (T) (curves 5 to in Fig. ||). 
Here it is reasonable to expect that V CI (T = T CI ) = and 
V cr grows monotonically to its zero temperature value as 
T — > 0. At T = (curve 0), R follows the power-law 
V~ as for V <C V CI (0). Figure H summarizes the second 
main result of the paper. 

Discussion: The present analysis may be tested with 
single- wall CNTs. Although a single- wall CNT has two 
transport channels rather than one, it can be analyzed in 
a similar way by exploiting the mapping from a single- 
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wall armchair CNT to a two-leg ladder |L6|. The con- 
ductance calculations closely parallel those for the sin- 
gle channel system. For Ug = 0, G crosses over from 
G NL = 4*T C+ /(1 + K C+ )G (T > T c ) to G SL = AK C+ G 
(T < T c ) Q. For small Uq, the power-law correction 
to G appears, whose temperature scaling exponents 0:1,2 
have new values, 01 = (K c+ — 1)/2(K C+ + 1) < and 
a 2 = K c+ — 1 < 0. For (effectively) large Uq, the lead- 
ing contribution to G again scales with temperature but 
the exponents 0:3,4 are modified to 0:3 = — 1 > 
and a 4 = {K~+ — l)/4 > 0. Then by interpolating the 
behaviors in the high and low energy scales, we find that 
the results in Fig. [I] remain equally valid for the two- leg 
ladder system. It can also be verified that the dV/dl-V 
characteristic remains qualitatively similar. 

Two recent experiments [T^J reported measurements 
on single- wall CNTs in contact with two superconducting 
leads. Thus the experimental configuration amounts to a 
S-CNT-S junction rather than a S-LL junction, and the 
Josephson current is observed in Ref. ||. In Rcf. [[To|| , 
however, the Josephson current is not observed and it is 
suggested that the system for some unknown reason be- 
haves as an "incoherent" series of two S-CNT junctions, 
where single S-CNT junction dynamics essentially deter- 
mines properties of the whole system. Motivated by this 
observation, we compare our main results (Figs.0, and|) 
with the results in Ref. ]ic|j . For a sufficiently clean junc- 
tion interface, the experiment found that the zero bias 
resistance shows nonmonotonic behavior, dR/dT > for 
T cr (« 4 K) < T < T c (w 10 K) and dR/dT < for 
T < T cr . This behavior is in qualitative agreement with 
the curves B and C in Fig. [l]. The experiment also found 
that the dV/dl-V characteristics above and below T cr 
are qualitatively different, again in qualitative agreement 
with Fig. |. 

However it is yet premature to interpret Ref. |l0[ ] as 
an experimental verification of the predictions for a S- 
LL junction. As recently pointed out by Wei et al. 
the level spacing of CNTs due to the finite size effect in 
Ref. jio) is of the same order as A. The level spacing 
can serve as a new energy cut-off and prevent the sys- 
tem from reaching the low energy regime in Figs, [l] and 
||. Moreover it is demonstrated Jl7|, [l^] that a noninter- 
acting model of a S-CNT-N junction, which takes into 
account the finite size effect, can produce nonmonotonic 
dependences. In this explanation, however, the resistance 
does not follow the power-law even sufficiently below T cr , 
in clear contrast to the predictions for a S-LL junction. 
Thus for an unambiguous verification of the predictions 
in this paper, we suggest a search for power-law behaviors 
in new experiments with longer CNTs, where the finite 
size effect is weaker. 

Summary: It is demonstrated that the resistance of 
a superconductor-Luttinger liquid junction with a suffi- 



ciently clean interface shows nonmonotonic temperature- 
and voltage-dependence due to the competition between 
the superconductivity and the repulsive interaction. The 
results are discussed in connection with a recent experi- 
ment ]l0| on single- wall carbon nanotubes in contact with 
superconductors. 
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